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ABSTRACT 

Context. Previous numerical studies have identified "phase mixing" of low-frequency Alfven waves as a mean of parallel electric 
field amplification and acceleration of electrons in a coUisionless plasma. 

Aims. Theoretical explanations are given of how this produces an amplification of the parallel electric field, and as a consequence, 
also leads to enhanced coUisionless damping of the wave by energy transfer to the electrons. 

Methods. Our results are based on the properties of the Alfven waves in a warm plasma which are obtained from drift-kinetic theory, 
in particular, the rate of their electron Landau damping. 

Results. Phase mixing in a coUisionless low-/? plasma proceeds in a manner very similar to the visco-resistive case, except for the 
fact that electron Landau damping is the primary energy dissipation channel. The time and length scales involved are evaluated. We 
also focus on the evolution of the parallel electric field and calculate its maximum value in the course of its amplification 
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1. Introduction 

At finite wave-numbers in the direction perpendicular to the am- 
bient magnetic field, Alfven waves produce a compression of the 
plasma which results in the creation of a parallel electric field 
via the thermo-electric effect, i.e. due to electron pressure fluc- 
tuations along the magnetic field lines. This electric field, whose 
magnitude increases with kj^, leads to wave -particle interactions, 
and hence, to coUisionless damping of the wave. 

Importance of this parallel electric field was pointed out 
already some time ago by (Hasegawa & Chen 1976). Indeed, 
they argue that "resonant absorption" (Hasegawa & Chen 1974) 
is a manifestation of mode conversion from the MHD Alfven 
wave (AW) to the kinetic Alfven wave (KAW) and that the 
physical mechanism of the heating depends on the collision- 
less absorption of the KAW. Although the original motivation 
was heating electrons in laboratory fusion plasmas, this electric 
field was also proposed as a mechanism which ca n accelerate 
electrons in space plasmasjH asegawa 1976; Haseg awa & Mimal 
1197 8'; 'Hasegawa| [T985l; lOoertz & Boswell 1979D and for under- 
standing solar coronal heatingj lonson 197^! 

(He vvaerts & Priest Il983h also introduced the idea of 
"phase-mixing" to improve the efficiency of AW dissipation. 
Their theory is based on visco-resistive magnetohydrodynamics 
(MHD). Since then, MHD phase mixing has attracted a signifi- 
cant amount of attention in the c ontext of heat ing open m agnetic 
struct u res in the solar coronal Parked 19911; TN akariak ov et al 
19971 iBothaetalJ l2000t iDe Moortel et'aiTl2000i; iHood et al 
2002) ■ Popular excitation mechanisms for coronal AWs in 



open magnetic structures are photospheric motions and chromo- 
spheric reconnection events, respectively for the low-frequency 
and high-frequency range of the spectrum. 

Phase mixing can be understood as the refraction of the 
wave while it propagates along a magnetic field with trans- 
verse variation in the Alfven velocity, i.e. the progressive in- 



crease of its kj _. This is a special occ urrence of conservative 
energy cascade dBian & Tsiklaurill2008l) . a phenomenon gener- 
ally attributed to non-linear interactions. Therefore, it is not sur- 
prising that phase mixing produces amplification of the parallel 
electric field that accompanies the Alfven wave in a coUisionless 
plasma, although this cannot be understood within the frame- 
work of ideal MHD theory which assumes E\\ - 0. 

Previous numerical studies of phase mixing in a coUisionless 
plasma have identified its implication in the generation of a par- 
allel el ec tric field and acceleration of electrons ( Tsiklauri et akj 
l20 05allbl:lTsikIauri & Harukill2008h . see also (Genot et aL, 19991 
2004) in the magnetospheric context. As stated above, the same 
features were established already some time ago, by Hasegawa 
and Chen, for resonant absorption. Here, we provide a detailed 
discussion of the role played by phase mixing in both parallel 
electric field amplification and enhanced electron Landau damp- 
ing of AWs in a coUisionless plasma. 

The calculations are based on the drift-kinetic theory pre- 
sented in Section II, which is valid in the limit of low-frequency 
fluctuations with a> «: cOci, cJd being the ion cyclotron frequency. 
Phase mixing and enhanced electron Landau damping of AWs 
in a coUisionless low-jS plasma are considered in Section III. 
Parallel electric field amplification is analyzed in Section IV. 
Conclusions and discussions are provided in in Section V. 

2. Kinetic properties of tlie Alfven wave in a warm 
coUisionless plasma 

Our starting point is the linearized drift-kinetic equation for the 
electrons: 



dtfi + V||V||/i - — £||5„„/() = 



(1) 



The latter is supplemented by Maxwell's equations. Faraday's 
law is 
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(j) is the electric potential. Ay the parallel component of the vector 
potential. The parallel component of Ampere's law reads 



ViA| 



Ane r 
= — J V,,/, 



\dv\\. 



(3) 



The above system is closed by the quasi-neutrality condition, 
which in the limit k^pi «: 1 , reads 



I 



fidv\ 



(4) 



Pi is the thermal ion Larmor radius at the temperature To, and 
«o is the background density. We set the Boltzmann constant to 
unity which means that the temperature has the unit of energy. 
While this so-called gyrokinetic Poisson equation [Eq.(4)] in- 
cludes the effect associated with the perpendicular ion polariza- 
tion drift, the electron response along the perturbed field lines is 
described by the drift-kinetic equation [Eq.(l)]. 

We assume a small deviation f\ from an equilibrium 
Maxwellian distribution /q: 



/o(v||) = 



no 



(5) 



The electron density perturbation is defined as rie — J /it/vy and 

the parallel current perturbation as J\\ - -e J V||/i(iv|| = -e«()M||e, 
u\\e being the electron parallel velocity, and the electron pressure 
perturbation is defined as Pg - nie J v^^fidv^ Hence, Ampere's 

law and Poisson law can be written respectively as V^^Ay - 
—(Alnc)J\\ and p^jV\e(j)IT(,i - ne/riQ. On one hand, taking the 
zeroth order moment of the electron kinetic equation provides 
the electron continuity equation: 



drie ., 

— + «OiK||M||, 







(6) 



On the other hand, the first moment provides the parallel electron 
momentum equation : 



name 



-ik\\Pe - noe-Eii, 



(7) 



It is usual to refer to the last equation as the Ohms's law and 
Pg = neTi)e for an isothermal plasma. Therefore, there are two 
possible sources of parallel electric field associated with the 
electron dynamics : inertia and pressure (or density) variations 
along the field lines. The continuity equation combined with 
Poisson law, yields a vorticity equation : 



d c(b c 

— p V — ^ H 



(8) 



Neglecting first the effects of electron inertia and electron pres- 
sure gradient in Ohm's law yields the MHD Ohm's law = 0, 
i.e. 



c dt 



(9) 



Introducing the stream and flux function for the velocity Uj. — 
xxVj_(p, and the magnetic field B^/ yfAjmomi = zx'Vj_4r, defined 
asifi = (c/B{))(p and i// = -Ay/ ^lAnnnii, gives 



at 



(10) 



(11) 



with va - Bq/ V47r«o'«i being the Alfven velocity. These two 
equations are the standard linearized reduced-MHD equations 
describing shear- Alfven waves with frequency ; 



(12) 



In the case where the parallel electric field is produced by den- 
sity fluctuation in Ohm's law, we have E\\ = -ik\\Toeine/no). 
Using the Poisson equation, E\\ = -ik\ip^^V\(p, which also re- 
veals the vortical nature of the parallel electric field. The pa- 
rameter ps = c jojci - yjTaJTQiPi is the ion gyroradius at the 
electron temperature. By including this parallel electric field in 
Ohm's law, an extension of the previous reduced-MHD system 
now takes the form 



d 2 

dtfr 2 2 

— = va/^||(^P-P,V^(/j), 



(13) 



(14) 



which describes the dynamics of kinetic Alfven waves with fre- 
quency 



(15) 



It is worth noticing that equations (fT3])-(fT4li can also be obtained 
directly from two-fluid MHD theory by retaining the Hall and 
electron pressure effects in Ohm's law(Bian & Tsiklauri 200^. 
Using the above results, it is easily seen that for kinetic Alfven 
waves, the magnitude of the parallel electric field is related to 
by 



17 ^^1 

En - — k\\ 



(16) 



IPs 



The above fluid derivation of the Alfven wave frequency 
gives the same result as its kinetic counterpart, however the lat- 
ter, which is presented below, is more complete in the sense 
that it also provides the imaginary part associated with Landau 
damping. The electron kinetic equation can be solved for the 
perturbed distribution function /i, i.e. 



(17) 



/i — ' — ^11 ^ < 

Me V7r^||V;^ V\\/Vte - (^lk\\V,e 

Some notations are introduced : x - v\\lvte, a - a>/k\\v,e and 
I r x" 2 

Z„(a) = — e-' dx, (18) 

y/n J X- a 

with Zo(a) being the standard plasma dispersion function. We 
also summarize some properties of the functions Z„\Z\ - 1 H- 
QfZo, Z2 = aZ\ . Moreover, in the limit o- <K 1 



Zo(a) ~ -la + i^{\ - a^). 



(19) 



Using the above properties, it follows that the density and current 
perturbations are related to the parallel electric field through: 



hdv\\ = [1 H- aZo(a')]£'||, 



for the density, and 



2/enow 

hn^n = — TTtCI -I- aZo(a)]£'||. 



(20) 



(21) 
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Hence, the relation between parallel current and parallel electric 
field is 



-10) 

4^ 



[1 + QrZo(a)]£'||. 



(22) 



It is convenient to define a coUisionless plasma conductivity a 
as 



7ii = (tE\\ 



(23) 



Its imaginary part results in the dispersion of the Alfven wave 
and the its real part yields the colhsionless dissipation. In the 
limit a = o)/k\\Vte <K 1, the real part is 



9 1/2 7 

e riooj 

"ll ^ Oe 



(24) 



This also gives the energy per united time transferred to the elec- 
trons through the relation : 



Q = ReiJiiEl) 



I.e. 



Ue 



(25) 



(26) 



with Aoe being the electron Debye length and Ue^^ =\ E^^ \ l%n 
being the energy density of the parallel component of the electric 
field. It is in fact a standard result that the asymptotic » 1 
averaged power transferred to electrons, Q - / vy < -eE\\f\ > 
dv\\ due to the presence of a an harmonic electric field fluctuation 
E\\ = cos(A:||Z - cot) is 

Q = -^:r^[v\\-^\=./k,. (27) 

This can easily be verified from Eq.(l) and for a Maxwellian dis- 
tribution it is equivalent to Eq.(26). Using the relation between 
E\\ and Bj_, Q can finally be expressed in term of the magnetic 
energy, Ub =\ \ /Stt, 

^ ^(J- k]_p] 



(28) 



The coefficient of proportionality between Q and U b, which has 
the dimension of the inverse of a time, is nothing else than the 
Landau damping rate. 

The Landau damping rate is now directly obtained, without 
any reference to its physical meaning, from the complex disper- 
sion relation. The kinetic dispersion relation is obtained from : 
ViA|| = jw/(A|4^c)[l+aZo(Qr)]£||,p?Vi0 = +aZo(a)]£|| 
and E\\ = -ik\\</> + io)A\\/c. It is 



piki + (1 - -^)[1 + aZoia)] = 0. 



VA 



(29) 



This is the general complex dispersion relation for the dispersive 
Alfven wave. In the limit a ^ 1, it reads 

a/ = k^vlil + kiplil - i ^fna)]. (30) 

Its real part corresponds to the frequency of the kinetic Alfven 
wave which was also derived from fluid theory above. Its imag- 
inary part, which corresponds to the Landau damping rate [see 
also Eq.(28)] and reads 



2 „2 



Most calculations above were finalized in the limit or <c 1 , in the 
opposite limit of a » 1 one obtains the frequency and damp- 
ing rate of the inertial Alfven wave, which has its parallel elec- 
tric field balanced by the electron inertia in Ohm's law. For fre- 
quency (jj ~ k\\VA, a ~ VAlvte, hence the kinetic Alfven wave 
regime corresponds to vaIvic 1 and the inertial Alfven wave 
regime to VaIvic » 1. In the foUowing we continue to focus on 
the warm plasma regime corresponding to 1 » j6e » me/rrii. 



3. Phase Mixing 

Phase mixing of a shear Alfven wave packet can be considered 
in the framework of an eikonal description: 



dx „ 

dk 

dt 



(32) 



(33) 



with 0) = ±k\\VA. These are the characteristics of the wave-kinetic 
equation 



^ -I- Vka>.V^e± - V^(oVke± = -y{k)e±. 
at 



(34) 



In the latter equation e± are the amplitudes of the wave-packets 
corresponding to w = ±k\\VA and -y(k) is a wave-number depen- 
dent damping rate. Following the trajectory of a wave-packet in 
phase space (x, k), its amplitude evolves according to : 



The latter equation is integrated to give 



e±(0 = e±(0)exp(- y(k)df). 



(35) 



(36) 



The principle of phase-mixing is simple: for any damping rate 
J which is an increasing function of k, any mechanism pro- 
ducing an increase in ^ as a function of time results also in a 
smaUer damping time scale. This is precisely the situation when 
the Alfven wave packet propagates along field lines with a trans- 
verse variation of the Alfven speed: the wave packet is sheared. 
In this case, say va(x) = -v^xz, z being the unit vector in the 
parallel direction and x the transverse coordinate, then 



dk± 
dt 



■ h^A^ 



(37) 



with by definition = vaIL^, being the characteristic length 
of the transverse inhomogeneity and = k\\{t - 0). This means 
that ki_ increases linearly with time due to diff'erential advection 
of the wave packets along the field fines, i.e. 



kAt) = kv'^t, 



(38) 



where we have taken kj_(t = 0) = without loss of generality. 

For a resistive MHD Ohms' law, E^ - rjJw the following 
results are well known. The damping rate is y(k) = ?7c(A^ -I- 
fcy)/47r - D,„{k\ + fcy), this is the fourier transform of the opera- 
tor responsible for magnetic diff'usion in the induction equation. 
Hence 



(31) e±(f) = e±(0)exp[-M I (1 + v'lt^)dt]. 



(39) 
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which in the limit f » ' yields 
e±(z)~e±(0)exp( 



(40) 



4. Parallel electric field generation 

For an Alfven wave created by a source through perturbation of 
the background magnetic field, a parallel electric field is pro- 
duced, provided kj^ is finite, which is given by Eq.(16); 



Since, z -VAt, we also have 



DmV'}u? , 

e±(0~ei(0)exp( -^z\ 



va , 



(41) 



for an Alfven wave excited at z = with frequency cj. In a col- 
lisionless plasma, when the dissipation is provided by electron 
Landau damping, with damping rate y(k) - ■\/7Tv\k\\k]_p'^/2v,e, 
the equivalent expressions are: 



e±(t) = e±(0)exp(-^^p2^^f3), 

6 V,e ■ " 



and also 



(42) 



(43) 



for an Alfven wave excited at z = with frequency a>. Hence, 
the phase mixing time scale is 



"^pm ~ 2/3, ' 

and the phase mixing length scale is 



(44) 



Ij-i 1/3,2/3 
V ^te W 

2/3 
Ps W 



(45) 



Notice that the scaling of the phase mixing length scale with the 
frequency oj in the spatial problem is different from that of re- 
sistive MHD phase mixing since the collisionless conductivity 
associated with electron Landau damping depends on a>, con- 
trary to Spitzer conductivity. However, the dependence with time 
or distance of the decay law, like exp(-Q'if^) or exp(-a'2Z^) are 
similar to resistive MHD phase mixing. The physical reason is 
obviously the common scaling of the damping rate y(k) with kj^ 
in the collisional and collisionless case. 

The following comments are due. The enhanced elec- 
tron Landau dampin g associated with phase-mixing was 



^pha 

first considered by (.Voitenko & GoossensI l2000a ). They de 
rived a relati on identical to Eg. (45) [se e equations (30) 
and (11) in dVoitenko & GoossensI l2000a')1. Moreover, re- 
sults of the Particles-In -Cell (PIC) simulations carried by 
(iTsiklauri & H aruki 2008') have produced /p„, oc oj^^ with ^ ^ 
1.10, for the dependence of the phase mixing length scale Zp,„ 
with frequency co. They also report that the parallel electric field 
associated with the Alfven wave is primarily balanced by the 
electron pressure gradient in their simulations. They attribute 
the scaling of /p,„ with oj to the effect of an "anomalous resis- 
tivity" due to "scattering of particles by magnetic fields" which 
"plays an effective role of collisions". Here, we emphasize that 
their PIC simulation results can be accurately interpreted as the 
"normal" effect of electron Landau damping of the KAW since 
it gives lp„, oc ai^^ with ^ = 1 . We now elaborate on the parallel 
electric field amplification which is observed in the simulations. 



(46) 



It is this parallel electric field which is responsible for the 
Landau damping of the wave (see above). For a given k\\ and 
6B^, this electric field is amplified provided that the kj_ asso- 
ciated with the wave field is also amplified, being a mono- 
tonic increasing function of kj_. However, E\\{kj_) also reaches a 
plateau for k^ps ~ 1, which is the boundary between the MHD 
and the dispersive regime. Indeed, 

£|l = ^-^k^k^plB^ (47) 

for kj^ps <K 1 and E\\ reaches its maximum, of the order of 

£|l = — ^iip.Bx (48) 
c 

when kj_ps ~ 1 or larger Therefore, significant amplification of 
this parallel electric field can only occur in the range of wave- 
numbers where the wave is non-dispersive, i.e. it behaves as a 
shear- Alfven wave with frequency a> - +k\\VA and, hence, it can 
be subject to standard phase mixing. 

From the results of the previous section we obtain the depen- 
dence with time of the parallel electric field strength during the 
phase mixing process : 



£||(f) = v'^k.pjexpi 



2/2 

A^A 2,3,3n 



(49) 



where a normalized electric field E\\ - E\\/(Bj_(0)va/c) has 
been defined. The variation with time E{t) has the form 
exp(-Q'if^), with a growth phase followed by a decay phase 
typical of the alternating field aligned current during phase mix- 
ing. Since z - vaI, then 



,2 ^2, 



E\\(z) 



■ exp(- 



VAVte 



-pWz\ 



(50) 



which has the form y82zexp(-Q'2Z^), for an Alfven wave excited 
at z = with frequency co. The above defined phase mixing 
time/length scales are precisely the scales associated with the 
amplification of the parallel electric field, i.e. the time/length 
scales for the parallel electric field to reach its maximum value 
given by : 

1/3 4/3 
~ (jJV.L pJ 

with 0) - k\\VA 



5. Conclusions 

Previous PIC (Particles-In-Cell) simulat ions of "colhsonle ss 
phase mixing" of Alfv en waves (Ts iklauri et alJ l2005allbL 
ITsiklauri & H aruki'2008') have identified its relation to the gen- 
eration of a parallel electric field and acceleration of electrons. 
Importance of this parallel electric field was first pointed out by 
Hasegawa and Chen in the context of resonant absorption, who 
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also showed that the dominant energy dissipation of the Alfven 
wave, in a collisionless low-/? plasma, involves energy transfer 
to the electrons(Hasegawa & Chen 1976). The role of electron 
Landau dam ping in "co llisonless phase mixing" was also con- 
sidered by ( V oitenko & Goossensi.2000a.b) 

Focusing on the "kinetic" regime of the dispersive Alfven 
wave, when v^/v/g <K 1, we provided a detailed discussion of the 
role played by phase mixing in both parallel electric field ampli- 
fication and enhanced electron Landau damping of the wave. 

Qualitatively, the physics of collisionless phase mixing can 
be summarized as follow. A parallel electric field accompa- 
nies the propagation of Alfven waves with finite k± because 
they compress the plasma. The magnitude of this electric field 
is an increasing function of k± that saturates in the dispersive 
range when kj_ps ~ 1 or larger. Therefore, any mechanism 
that produces an increase in kj_ also leads to the amplification 
of the parallel electric field associated with the Alfven wave. 
Phase mixing is such a mechanism, independently of the en- 
ergy dissipation channel. Phas e mixing is a special oc currence 
of energy-conserving cascadel' Bian & Tsiklaurill2008h . Such a 
cascade, predominantly involving perpendicular wave-numbers, 
is generally attributed to non-linear interactions, i.e. to turbu- 
lence. Existence of this parallel electric field and the depen- 
dence of its magnitude with yield a Landau damping rate 
which scales like k]_, just as visco-resistive damping. This can be 
demonstrated very simply in the framework of drift-kinetic the- 
ory. Therefore, in a collisionless plasma, phase mixing leads to 
enhanced electron Landau damping of the Alfven wave in a man- 
ner which is very similar to the well-studied case of enhanced 
visco-resistive damping. As a consequence, once the wave has 
Landau damped in a collisionless low-^S plasma, its energy has 
been transferred to the electrons, and the time and length scales 
involved have been evaluated for small amplitude perturbations. 
Moreover, we studied the evolution of the magnitude of the par- 
allel electric field in the course of its amplification and calculated 
its maximum value. 

We argued that the scaling of the phase mixing length 
scale with fr equency, j™ oc a»"^ and ( ^ 1, reported by 
dTsiklauri & H aruki 2008|) has a simple interpretation in term 
of electron Landau damping. PIC simulations of collisionless 
phase mixing are valuable tools because they can provide direct 
information on the modification of the electron di stribution func- 
tion involved in the acceleration process, see dTsiklauri et al.l 
[2005a.b), a feature which the present kind of analysis is not ca- 
pable of. For this, a theoretical framework is needed, e.g. quasi- 
linear theory. This is the subject of ongoing work. 
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